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ABSTRACT: In this study, we developed kinetic models to predict the growth of pathogenic Escherichia coli on cheeses during 
storage at constant and changing temperatures. A five-strain mixture of pathogenic E. coli was inoculated onto natural cheeses (Brie and 
Camembert) and processed cheeses (sliced Mozzarella and sliced Cheddar) at 3 to 4 log CFU/g. The inoculated cheeses were stored at 4, 
10, 15, 25, and 30°C for 1 to 320 h, with a different storage time being used for each temperature. Total bacteria and E. coli cells were 
enumerated on tryptic soy agar and MacConkey sorbitol agar, respectively. E. coli growth data were fitted to the Baranyi model to 
calculate the maximum specific growth rate (/^max; log CFU/g/h), lag phase duration (LPD; h), lower asymptote (log CFU/g), and upper 
asymptote (log CFU/g). The kinetic parameters were then analyzed as a function of storage temperature, using the square root model, 
polynomial equation, and linear equation. A dynamic model was also developed for varying temperature. The model performance was 
evaluated against observed data, and the root mean square error (RMSE) was calculated. At 4°C, E. coli cell growth was not observed on 
any cheese. However, E. coli growth was observed at 10°C to 30°C with a ^^^^ of 0.01 to 1.03 log CFU/g/h, depending on the cheese. 
The /i^ax values increased as temperature increased, while LPD values decreased, and /(,„,x and LPD values were different among the 
four types of cheese. The developed models showed adequate performance (RMSE = 0.176-0.337), indicating that these models should 
be useful for describing the growth kinetics of E. coli on various cheeses. (Key Words: Escherichia coli. Cheese, Predictive Model, 
Dynamic Model) 



INTRODUCTION 

Cheese is one of the most popular foods globally. The 
International Dairy Foods Association (IDFA) reported in 
2010, that cheese is the major manufactured dairy product, 
and that its importance for the dairy industry has grown 
remarkably. The consumption of cheese has been increasing 
gradually since the 1990s (Kim et al., 2007). However, 
some cheese is occasionally contaminated with food-borne 
pathogens such as Listeria monocytogenes. Staphylococcus 
aureus, and Escherichia coli (E. coli, Thayer et al., 1998; 
Kaan Tekin§en and Ozdemir, 2006; Jo et al., 2007). 

E. coli, a facultative anaerobic gram-negative bacillus, 
is commonly found in the intestinal flora of man and 
animals, and certain strains are pathogenic (Olsvik et al.. 
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1991). Pathogenic E. coli strains are classified by infection 
and pathologic mechanism: enteropathogenic E. coli 
(EPEC), enteroinvasive E. coli, enterotoxigenic E. coli 
(EXEC), enterohemorrhagic E. coli (EHEC), and 
enteroaggregative E. coli (EAEC) (Nataro and Kaper, 1998). 
Recently, there have been reports from many countries, 
regarding the isolation of E. coli from various cheeses 
(Haran et al., 2012; Zinke et al., 2012). Moreover, several 
outbreaks of E. coli infections, related to cheese, have been 
reported (FSN, 2010). Therefore, several countries have a 
quantitative standard or "zero tolerance" policy for 
controlUng pathogens in cheese (Health Canada, 2008; 
FDA, 2009). 

Predictive models can be used to determine how 
intrinsic or extrinsic factors affect or interact with growth 
parameters (Lihono et al., 2003). Moreover, they can 
predict the growth or survival of bacteria in accordance 
with a diverse range of parameters, such as storage 
temperature and storage time (McMeekin et al., 1997). 
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Storage temperature is a major factor for bacterial growth, 
and it can be changed continuously under various 
circumstances such as transportation, display at retail, and 
storage at home. Therefore, it is important to describe the 
growth kinetics of pathogens under dynamic storage 
conditions. 

Although the growth of E. coli has been evaluated in 
raw milk and unpasteurized cheese using predictive models 
(Sutherland et al., 1995), the growth of pathogenic E. coli 
on cheese has not been studied, especially in response to 
varying temperature. Therefore, the objective of this study 
was to develop a mathematical model that describes the 
growth kinetics of pathogenic E. coli under constant and 
dynamic temperature conditions. 

MATERIALS AND METHODS 

Inoculum preparation 

Five strains of E. coli, NCCP14037 (ETEC), NCCP 
14038 (EPEC), NCCP 14039 (EAEC), NCCP 15661 
(EPEC), and ATCC11142 (ETEC) were suppHed from 
National Culture Collection for Pathogens. They were 
cultured in 10 mL of tryptic soy broth (TSB; Bacto, Becton 
Dickinson, MD, USA), and incubated at 35°C for 2 h. We 
inoculated 0.1 mL of each culture suspension into 10 mL of 
TSB for subculturing at 35 °C for 24 h. After incubation, the 
subcultures were then mixed and harvested by 
centrifugation for 15 min at 1,912 g and 4°C. The cell 
pellets were thoroughly washed twice with phosphate- 
buffered saline (PBS, pH 7.4; 0.2 g/L KH2PO4, 1.5 g/L 
Na2HP04, 8.0 g/L NaCl, 0.2 g/L KCl in distilled water), 
followed by serial dilution with PBS to obtain 
approximately 5 to 6 log CFU/mL. 

Sample preparation and inoculation 

Two commercial natmal cheeses (Brie and Camembert; 
no antimicrobial included) and two processed cheeses 
(sliced Mozzarella and sliced Cheddar; no antimicrobial 
included) were purchased at a retail store. The natural 
cheeses were cut into 15 g portions, and transferred into 
sample bags. The surfaces of the cheeses were inoculated 
with 0. 1 mL of inoculum at 3 to 4 log CFU/g. The cheese 
samples were then massaged to spread the bacteria, and 
sealed (Food Guard VP5700, RoUpack, Gyeonggi, Korea). 
In addition, 0.1 mL portions of inoculum were inoculated 
onto each slice (18 g) of processed cheese (Mozzarella and 
Cheddar), using a sterile bent spreader. The processed 
cheese samples were then covered aboriginally. Two 
samples were then placed in a plastic bag, and sealed, 
followed by storage at 4°C (1,320 h), 10°C (768 h), 15°C 
(120 h), 25°C (48 h), and 30°C (24 h). 



Experimental growth analysis 

Cheese samples were analyzed 9 to 16 times during 
storage, depending on the storage temperature. Thirty 
milhliter of 0.1% buffered peptone water (BPW; Difco, 
Becton Dickinson, MD, USA) was added to the natural 
cheese samples and the mixture was homogenized for 120 
min in a blender (BagMixer; Interscience, France), while 
processed cheese samples were transferred into a filter bag 
(Sample bag, 3M, Korea), containing 30 mL of BPW, and 
blended for 2 min. The homogenates were then serially 
diluted with BPW. Aliquots (0.1 mL) of the diluents were 
surface-plated on tryptic soy agar (TSA; Difco) and 
MacConkey sorbitol agar (Mac; Difco) for determining 
total bacterial and E. coli counts, respectively. The plates 
were incubated at 35°C for 24 h, and the colonies were 
manually counted. The pH values of the homogenates were 
measured via a digital pH meter (Accumet, Denver 
Instruments, CO, NY, USA), and the water activity, a„ of 
the samples was determined by a water activity meter 
(Aquaspector, NAGY Messsysteme, Gaufelden, Germany). 

Primary model 

The E. coli growth data was fitted to the Baranyi model 
(Baranyi and Robert, 1994), using DMFit (Institute of Food 
Research, Norwich, UK) to estimate the maximum specific 
growth rate (/^maxi log CFU/g/h), the lag phase duration 
(LPD; h), the lower asymptote (A^o; log CFU/g), and the 
upper asymptote (A^max; log CFU/g). 

Secondary model 

The kinetic parameters from the primary model were 
further analyzed as a function of storage temperature. The 
square root model and polynomial model were used for 
natural cheese and processed cheese, respectively, as 
follows: 

fl,x(r-T^„) (1) 

/^max = ai+a2X.T+aiX.'f (2) 

where is the slope of the linear regression for the 
square root of fia^mA Tis the temperature (°C). T^,, is the 
theoretical minimum temperature, a, is the coefficient of a 
polynomial equation. In addition, a linear equation was 
fitted to the square root of the inverse LPD for both natural 
and processed cheeses using the expression, 

= +«LPD X«LPD (3) 

where are the slopes of the regression lines for the 
square root of the inverse LPD, and Tis temperature (°C). 
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Validation 

To validate the models, E. coli cell counts were obtained 
experimentally. These data were then compared with the 
predicted E. coli cell counts, which were estimated by our 
models. Subsequently, the root mean square error (RMSE) 
was calculated to evaluate the model performance as 
follows: 



RMSE = 



1^ (observed values - predicted values)^ 



(4) 



where n is the number of observations. 

Dynamic model 

A study by Lee et al. (2008) reported that the mean 
temperature and storage time in a retail store were 7°C (Tj^m, 
-2°C; rn,„x, 22.9°C) and 2 days, respectively, and that 
purchased cheeses were then transported in a car at 18°C 
for up to 1 h. In addition, Bahk (2010) showed that the 
mean temperature and storage time for a home refrigerator 
were 4°C (T^,n, -5°C; r,nax, 14°C) and 10 days, respectively 
Therefore, E. coli populations were simulated by using the 
mathematical model defined by Baranyi and Roberts (1994) 
and a temperature profile based on the above studies (4°C 
to 14°C). 

Statistical analysis 

The experiments were repeated twice with two samples 
per repeat (n = 4). The growth parameters (/imax, LPD, Nq, 
and A'max) were analyzed using the general Uinear model 
procedure of SAS version 9.2 (SAS Institute, NC, USA). 



Mean comparisons were performed by using a pairwise t- 
test with the criterion for significance set at p<0.05. 

RESULTS AND DISCUSSION 

pH values were similar for the four types of cheese and 
bacterial growth was not affected by pH (data not shown). E. 
coli growth was not observed for any cheese type at 4°C 
(data not shown). The processed Cheddar cheese showed E. 
coli growth at 15°C to 30°C, and the other cheeses had 
bacterial growth at 10°C to 30°C. LPD values decreased as 
the storage temperature increased (p<0.05), but no 
significant difference was observed between the cheeses 
(Table 1). /Zmax values increased (p<0.05) when storage 
temperature increased, and the values of natural 

cheeses were higher (p<0.05) than those of the processed 
cheeses (Table 1). The a^. values of Brie, Camembert, 
Mozzarella slice, and Cheddar shce cheese were 0.991, 
0.987, 0.975, and 0.973, respectively, bus, these lower /Uraax 
values of processed cheeses may have been caused by a 
lower a„ (0.973 to 0.975), compared with those (0.987 to 
0.991) of the natural cheeses. E. coli growth slows as a„ 
approaches its minimum (0.950) (Aberoumand, 2010). The 
low a^^ of processed cheeses may be caused by evaporation 
during the heating process; NaCl concentrations were 
similar among the cheeses. In addition, the A^max values were 
lower (p<0.05) in processed cheese than in natural cheese 
(Table 1). The Revalues of the primary models ranged from 
0.918 to 0.998, indicating that the fit between the primary 
model and the E. coli growth data was appropriate (Table 1). 

To evaluate the effect of storage temperature on kinetic 



Table 1. Kinetic parameters (mean±standard error) of pathogenic Escherichia coli on natural and processed cheeses, calculated by the 
Baranyi equation (Baran)d and Roberts, 1994) 



Cheese 




Storage 
temperature (°C) 


LPD (h) 


(log CFU/g/h) 


A'o 

(log CFU/g) 


^ * max 

(log CFU/g) 




Natural 


Brie 


10 


20.13+9.60^^ 


0.03+0.00'^ 


3.6+0.1* 


8.6+0.3" 


0.987 


cheese 


cheese 


15 


9.61+2.47'' 


0.07±0.01'^'^ 


3.5+0.0* 


9.4±0.2* 


0.992 






25 


5.40+0.55'' 


0.4540.03" 


3.6+0.1* 


9.1+0.2*" 


0.973 






30 


5.68+0.46'' 


0.94+0.10* 


3.6±0.0* 


9.1±0.1*" 


0.996 




Camembert 


10 


24.49+3.85"^ 


0.03+0.00'^ 


3.4+0.1* 


8.4±0.3" 


0.998 




cheese 


15 


10.35±2.08'' 


0.09+0.01°'^ 


3.4+0.0* 


9.3+0.1* 


0.984 






25 


5.74±1.06'' 


0.44+0.03" 


3.4+0.0* 


9.1+0.1*" 


0.982 






30 


5.92+0.71" 


1.03±0.07* 


3.6+0.0* 


9.0+0.1*" 


0.985 


Processed 


Mozzarella 


10 


3.44+8.10* 


o.oi+o.oo'^ 


3.0+0.1" 


6.3+0.1™ 


0.918 


cheese 


slice cheese 


15 


9.22±4.45'' 


0.06±0.0l'"^ 


3.1+0.0" 


7.7+0.1™ 


0.984 






25 


4.30+1.88" 


0.25+0.05™ 


3.3+0.1" 


7.8±0.2° 


0.988 






30 


1.79+0.69" 


0.33+0.01^^ 


3.0±0.0" 


8.0±0.l'^ 


0.992 




Cheddar shce 


15 


35.95+16.22* 


0.03+0.00'^ 


2.8±0.0" 


7.2±0.6™ 


0.967 




cheese 


25 


6.83*1.15" 


0.18±0.01° 


2.9+0.2* 


7.340.0^^ 


0.977 






30 


6.()8±2.()4'-""' 


0.28±().()7'''' 


B.OiO.l" 


7.()±().l™ 


0.988 



LPD, lag phase duration; /ii„ax, maximum specific growth rate; A'o, lower asymptote, Nmax, upper asymptote. 
Different letters in a same column mean significantly different at p<0.05. 
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parameters such as /Vax and LPD, secondary models were 
developed (Figures 1 and 2). When there was no growth of 
E. coli, e.g., at 4°C, /ZmaxWas estimated to be zero, and the 
length of the storage period was estimated to be the LPD, 
which is the period of the time, showing no growth. The 
model predictions were close to the experimentally 
observed kinetic data, with high values (0.890 to 0.994), 
indicating that the secondary models were valid for 
describing the effect of storage temperature on kinetic 
parameters. 

To vahdate our model, predicted ^^ax and LPD values 
were calculated via the secondary model at specific 
temperatures (18°C, 23°C, and 28°C). These predicted 
kinetic parameters were then used to predict E. coli cell 
counts at given storage times in accordance with the 
primary model. The predicted E. coli cell counts were then 
compared with the observed E. coli cell counts, which were 
obtained from our experiments. The RMSE values of Brie, 
Camembert, Mozzarella, and Cheddar slice cheeses, 
indicating the distance between observed values and 
predicted values, were 0.264, 0.218, 0.176, and 0.337, 
respectively. , depending on the cheese. This result indicates 
that the model performance is adequate. The RMSE values 



for models by Perez-Rodriguez et al. (2013) and Skfivanova 
et al. (2008) were 0.30 to 1.33 and 0.29 to 0.78, respectively, 
and this was considered adequate by the authors. Lee et al. 
(2012) concluded that their model was adequate with 
RMSE values of 0.326 to 0.361. Therefore, the performance 
of the model described by this study can be considered 
adequate. 

E. coli growth on cheese was also simulated for 
changing storage temperatures, ho values, which indicate 
the initial physiological state of bacteria (Grijspeerdt and 
Vanrolleghem, 1999) were calculated to be 1.816 to 2.074 
for natural cheese and 0.549 to 0.874 for processed cheese 
(Figure 3). Although some experimental data were not very 
close to the predicted values, most data were within 95% 
confidence intervals. This result suggests that the dynamic 
models are useful for prediction of E. coli cell growth on 
natural and processed cheeses under dynamic temperature 
conditions. 

In conclusion, the models developed in this study can be 
used to describe E. coli growth kinetics on natural cheeses 
(Brie and Camembert) and processed cheeses (sliced 
Mozzarella and sliced Cheddar). Such models could prove 
useful for exposure analysis in microbial risk assessment. 
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Figure 1. Square root model for (A-B) and linear equation for lag phase duration (C-D) developed for Brie and Camembert cheeses. 
•, observed value; — , predicted Une; — , 95% confidence interval. 
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Figure 2. Polynomial equation for /Vax (A-B) and linear equation for lag phase duration (C-D) developed for Mozzarella slice and 
Cheddar shce cheeses. •, observed value; — , predicted line; — , 95% confidence interval. 
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Figure 3. Predicted Escherichia coli growth on Brie cheese (A), Camembert cheese (B), Mozzarella slice cheese (C), and Cheddar slice 
cheese (D) under dynamic temperature condition. •, observed data, — , predicted line, ,95% confidence interval, - • -, temperature. 
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